High-field limit of the Vlasov-Poisson-Fokker-Planck system: A comparison of 

different perturbation methods 



Luis L. Bonilla 

Escuela Politecnica Superior, Universidad Carlos III de Madrid, Avda. Universidad 30, 
28911 Leganes, Spain. E-mail: bonilla@ing.uc3m.es 

Juan Soler 

Departamento de Matemdtica Aplicada, Facultad de Ciencias, Universidad de Granada 
18071 Granada, Spain. E-mail: jsoler@ugr.es 
(June 30, 2000) 

A reduced drift-diffusion (Smoluchowski-Poisson) equation is found for the electric charge in the 
high-field limit of the Vlasov-Poisson-Fokker-Planck system, both in one and three dimensions. The 
corresponding electric field satisfies a Burgers equation. Three methods are compared in the one- 
dimensional case: Hilbert expansion, Chapman-Enskog procedure and closure of the hierarchy of 
equations for the moments of the probability density. Of these methods, only the Chapman-Enskog 
method is able to systematically yield reduced equations containing terms of different order. 



I. INTRODUCTION 

The Vlasov-Poisson-Fokker-Planck (VPFP) system de- 
scribes the probability density p of an ensemble of electri- 
cally charged Brownian particles in contact with a ther- 
mal bath and subject to their self-consistent electric field 
cj), see p[. It consists of the equations: 



dp 
dt 



pv H V^c 

v m 



pk B T 



V v p 



= 0. 



eA x (j> = -eN Q pdv 



p{x, v, t) dx dv = 1. 
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(2) 



(3) 



Here ,/Vo is the number of particles, p a friction coefficient, 
e and m particle charge and mass, T temperature and e 
permittivity. Sometimes it is interesting to consider the 
high- field limit, in which the collision term (divergence 
with respect to v) in (|l|) dominates the other two. In this 
case, p quickly evolves toward a displaced Maxwellian, 



which is the thermal velocity. If we impose that pv and 
eV ' x <fi/m be of the same order, we find the unit of electric 
field, Eq, 
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yj mk B T . 
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The Poisson equation yields the unit of length, bEq/xq 
eN Q /xl, or 



xq 



eE \pe/ 



(6) 



The relation between the second and third terms in (|l|) 
yields the small parameter 



c = 



px 



1 fe 2 (k B Tf 



e \m(pN()) 2 



(7) 



The unit of p should be (xqVq)~ 3 because of the normal- 
ization condition (^). Finally, the unit of time should be 
xq/vq provided we want the first two terms in (|l|) to be 
of the same order. Nondimensionalizing the system ([l]) - 
(ft) in these units, we can rewrite it as 



2nk B T 



exp 



m [v + — V 2 



2k B T 



times a slowly- varying function, P(x,t). Essentially, P 
is a charge density whose evolution we would like to 
describe. Appropriate scales in which this fast equi- 
libration to the displaced Maxwellian occurs may be 
found as follows. A typical velocity is found by bal- 
ancing friction, pvp and flux pk B TV v p/m, which yields 
pv a = pk B T/(mv ), or 



kuT 



rn 



(4) 



V, • + Y., <.•)/> + V,./>] = < (?£+vV x p) . !N) 



A x ct) = -Jpdv, (9) 
pdxdv = 1, (10) 



to be solved together with appropriate initial conditions 
and decay boundary conditions at infinity. 

We are interested in finding a simpler reduced equa- 
tion for P(x,t) in the high-field limit as e — > 0+. This 
problem has been recently tackled by Nieto, Poupaud 
and Soler, g], who proved rigorously that P obeys the 
following hyperbolic system (in dimension one): 



1 



p t = 1 ton 



*XX 

Pdx = l, 



(11) 
(12) 

(13) 



to leading order as e — > 0+ (subscripts mean partial 
derivation with respect to the indicated variable). By 
using the decay of P at infinity, they derived the follow- 
ing equation for the electric field, E = —(f) x : 



Et + E E T 



0. 



(14) 



As is well-known, this equation may develop shock waves 
in finite time. Nieto et al proved that there exists a 
unique entropy solution of the VPFP system in the high- 
field limit, and that the shock velocity corresponds to in- 
terpreting ( |l4| ) as the conservation law E t + ^(E 2 ) x = 0. 
The uniqueness result suggests that this conservation law 
could be regularized by adding a small viscosity term, 
e E xx . Then analysis of the resulting Burgers equation 
in the limit as e ^ 0+ yields the unique entropy solu- 
tion. Their method amounted to closing the hierarchy 
of equations for the moments of p, and leaves open the 
question of how to add higher-order terms to the system 
(ll]) - (p"3"|), or (which is related) how to regularize (|l4|). 

In this paper, we shall consider these problems. To 
this end, we shall compare three classical approaches in 
kinetic theory, the Hilbert expansion, closure of moment 
equations and the Chapman-Enskog method (CEM) as 
applied to the one-dimensional VPFP system. These ap- 
proaches yield the same result for the parabolic scaling 



V„ • (v p + V„p) =e(v- V x p - V x (f> ■ V v p) + e" 



dp 
dt 



(15) 



which is interesting and usually employed in kinetic the- 
ory . This scaling corresponds to lower values of the 
electric field, so that the diffusive and frictional terms 
in the VPFP equation are in fact much larger than the 
corresponding term in (Q). For the present model, the 
parabolic scaling was also studied rigorously by Poupaud 
and Soler p], who derived from ( |T5| ) a drift-diffusion for 
P = J pdv in the limit as £ ^ 0+. In the high- field 
limit corresponding to the scaling (^J), we shall see that 
only the CEM yields satisfactory results systematically. 
In fact, the Hilbert expansion yields the following addi- 
tional system for the order e correction to P, P' 1 ': 







p t-T x ^ )pil) 



(16) 
(17) 
(18) 



At face value, this equation breaks down in the shock 
regions. However, it is compatible with the result of the 
CEM, 



^P) = P^ 
•Si) = _p(l) 

'xx 1 1 

P^dx = 0. 



Pt - Q^{<l> x P)=eP xx , 



^XX " P-i 

Pdx = 1, 



(19) 
(20) 
(21) 



provided we substitute P by P + and equate like 

powers of e in both sides of these equations. This sit- 
uation is exactly that found when using the method of 
multiple scales to describe certain codimension two bi- 
furcations The key point is that the CEM yields 
reduced equations which may contain terms of different 
order in e, unlike the other methods. Appropriate closure 
of the equations for the moments of p may yield ( |l9| ) - 
(pl|), but it is not a systematic method. 

The rest of the paper is as follows. The Hilbert expan- 
sion and the CEM are applied to the one-dimensional 
VPFP system in Sections || and [II, respectively. The 
latter section also contains the result of applying the 
Chapman-Enskog procedure to the three-dimensional 
problem. Section [IV] shows how to obtain the same re- 
sults by closing the hierarchy of equations for the mo- 
ments. Finally Section ^ contains our conclusions. 



II. THE HILBERT EXPANSION 

It consists of inserting the power series 

p = p^+ep^+e 2 p^+0(e 3 ), (22) 

in the VPFP system. This expansion is really akin to the 
method of multiple scales with slow time scale t S. We 
assume that the leading order contribution to the solu- 
tion has already relaxed in the fast time scale. Equat- 
ing like powers of e in both sides of the one-dimensional 
VPFP system, we obtain 



dx 2 



p<® dv, 



p {0) dxdv = 1, 



(23) 
(24) 
(25) 
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(26) 
(27) 

p {1) dxdv = 0, (28) 
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poo 

W - - / dv, 



and so on. 

The solution of (|2 



j {2) dxdv = 0, 



3) is 



(29) 
(30) 

(31) 



P(z,i), 



for a function P to be determined and such that 

= -P, 



P(x,t) dx = 1. 



(32) 

(33) 
(34) 

(35) 



Notice that the Fokker-Planck collision term behaves 
more nicely than the general Boltzmann term. In the 
high-field limit, the latter may give rise to runaway solu- 
tions for many forms of the collision frequency Q . 

Since J Cp^dv = 0, ( p6| ) yields the following solvabil- 
ity condition 

f p (n) dv + JL [ vp Wdv = Q. 



Of 



dx 



necessary for p^> to be bounded even as v — > ±oo. This 
provides 



P. 



(36) 
(37) 



together with (|13|), which coincides with Nieto et aFs re- 
suit, §. 

Suppose we now want to correct the leading order re- 
sult by going one step further, to solving (26) - (p8|). The 
result is 



= (p«-[P x + ( ^)p] 



V 



V — 1 



p< 



(38) 

(39) 
(40) 



where P^(x,t) is yet to be determined. By using the 
solvability condition for (S), we obtain 



P (1) dx = 0, 



P 



(i) 



d_ 

dx 



(0)pW j. ^(i) P\ - 



P) = P X 



(41) 



to be solved together with (|39j) and (|40|). Notice that 
the left hand side of (fll|) is a linearization of (^6|) about 



P. This equation does not make sense at those points 
where P is discontinuous, but we can easily see that ( |3^ ) 
and ( pLl] ) are obtained by taking V ~ P + ePW in the 
following equation: 



Pt-^ {<t> x V + eP x ) =0(e 2 ), 



(42) 



to be solved together with the Poisson equation <fi xx = 
—V and J Vdx = 1 . This situation is exactly that found 
when using the method of multiple scales to find the nor- 
mal form describing certain codimension two bifurcations 
HQ . The Hilbert expansion is akin to the method of mul- 
tiple scales in that it yields reduced equations all whose 
terms are of the same order. Thus if we insist in calculat- 
ing higher-order equations, we obtain coupled systems of 
equations such as those written above. To find directly 
a reduced equation including terms of different order in 
e, we should use the CEM. Appropriate closure of the 
equations for the moments of p may yield (^2|), but not 
in a systematic manner. 



III. CHAPMAN-ENSKOG METHOD 



In one dimension, the dimensionless VPFP system is 



d 2 4> _ 

dx 2 



pdv, (44) 



to be solved together with 

hoo />+oo 



p(x, v, t) dx dv = 1, 



(45) 



oo J — oo 



and appropriate initial and decay conditions as v and x 
tend to ±oo. 



A. Chapman-Enskog method 

Setting e = in (|43"|), we find a simple equation to 
be solved together with ( p4| ) and (|45|). Its solution is a 
displaced Maxwellian: 



e 2 



POM), 



V27T 

V = v 

4°i = -P 

/oo 
P{x,t) dx = 1. 
-OO 



(46) 

(47) 
(48) 

(49) 



Notice that P(x, t) is an arbitrary function of x and 
t except for (M9). Furthermore, (Eq) correspond to a 



3 



particular choice of initial conditions. The Chapman- 
Enskog ansatz consists of assuming that p has the fol- 
lowing asymptotic expansion 



P(x,t;e) + 



n=l 



(50) 



Furthermore, we impose that the amplitude P obeys an 
equation: 



dP 
~dt 



(51) 



n=0 



where are functionals of P to be determined as the 
procedure goes on. This equation for P is not explicitly 
written in the usual presentations of CEM |||| . Instead, 
the form of this equation is guessed by writing equations 
for the moments of p and using gradient expansions. We 
find this latter procedure more confusing. 

Insertion of (|50| ) and (|5l]) into the equations and aux- 
iliary conditions yields a hierarchy of linear equations 
for the p( n > . Notice that the latter depend on time only 
through their dependence on P. The functionals F^ n '{P) 
are determined so that each equation (and set of auxil- 
iary conditions) for p^ has a solution which is bounded 
for all values of v, even as v — > ±00. Once a sufficient 
number of is determined, (|5l| ) is the sought am- 

plitude equation. Please notice that, unlike results from 
the method of multiple scales, terms in ( j5l|) may be of 
different order. 

Let us illustrate how the procedure works by finding 
F(°) and . Insertion of © and @ in (||), @ and 
( |l0| ) yields the following hierarchy of linear equations: 



Cp 



(i) - 
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(52) 



Cp 



(2) 
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P (1) - VP<j) 



(0)1 

xt \p t= F( 1 ) 
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(53) 



and so on. Before equating like powers of e, we have sub- 
stituted Pt = X)^o wnen performing time differen- 
tiations such as 4>ft- This yields the obvious terms in 
the hierarchy of equations, which results in rather clut- 
tered formulas, as we ascend in the hierarchy of equa- 
tions. These equations are to be supplemented by the 
normalization conditions (E^), 



p {n) dv = 0, n > 1, 



(54) 



Cp 



(«) = Ji 



d_ 

dv 



p 



/27T 



(55) 

(56) 



for n = 1, 2, V is again given by (p3|). 

Let us now consider (52). Since the v integral of its 
left hand side is zero, this equation has a solution only 
if the v integral of its right hand side is zero. The corre- 
sponding integrals are simplified by using the symmetry 
of the Maxwellian and shifting integration variables from 
v to V. The condition that the integral of the right side 
vanish yields 



P 



(57) 



Notice that we needed F^ ^ in the right side of (52) for 
this equation to have an appropriate solution. 

We now calculate 4> x °t in order to simplify the right 
side of (^2|). As we explained above, at this order we 
jp(o) - fjf» 



should substitute Pt 
Hence 



K >p }, 



when needed. 



= -K x *P t = -K x * 

= -K xx 



i°) p) = p 



(58) 



where K is the Green's function of the one-dimensional 
Laplacian and * means convolution product. In the pre- 
vious calculation, we have integrated by parts and used 
the decay properties of P and the symmetry of K. In- 
serting now ([57]) and ( j58l) in Eq. (p2]), we realize that 
the right hand side thereof is the partial derivative of 

(v^p-p x ) e - v2 /y 

mediately yields 



2tt with respect to v. This im- 



2?r 



VP a 



(59) 



which satisfies ( |54|) and yields (p^ — 0. In (|59|) , we have 
omitted adding a term like that proportional to PW in 
the right hand side of (|3|) (satisfying J P^ 1 ) dx — 0) be- 
cause all such terms are already included in the amplitude 
P(x,t;e). 

To find F^\ we insert ( |59| ) in ( p3| ) and use the solv- 
ability condition for this equation. Simplifications arise 
from the identities 



p {1) dv = 0, 



and 



VPx' 



dx 



vp 



« dv. 



and the linear equations and definitions: 



The result is 



4 



= P, 



(60) 



The solution is then 



We can now insert (|57j ) and (|6C|) into (|5lJ) to obtain 
the sought Smoluchowski equation for P: 



P {1) = 



(61) 



(27T) 



to be solved together with (|4S|). Substituting ( ^g ) into 
(|6l"|), and integrating with respect to x, we find the fol- 
lowing Burgers equation for E = — <p x : 



~(VV-I): V X V X 



PX7 X ; 



(69) 



Et + E E x — e E x 



(62) 



Provided that the initial electric field be uniformly 
bounded, this Burgers equation has solutions which con- 
verge to the distributional solutions of the Hopf equa- 
tion (|lj). Thus the CEM is in fact a vanishing- viscosity 
method which yields a unique solution approaching the r 
unique entropy solution of @ as e -> 0+. In Ref. |§, P (1) = V x ■ I PV X 
the authors proved convergence towards the unique en- ^ 
tropy solution of the Hopf equation by means of a differ- 
ent method which used that the limiting electric field is 
a monotone decreasing function. 



Notice that in the one-dimensional case, the term 
V x [cj>f^ — \ v *^ - J vanishes. We could have added a 
term which solves the homogeneous equation Cu = in 
the right hand side of (69). However, as we explained 
before, the effect of this term is already included in 
P(x, t; e). The solvability condition for the p^ equation 
yields 



,(0) 



V r P 



(70) 



The reduced Smoluchowski equation for P is therefore 



Pt = V x -iPV x <f> + ePV x 



B. Three-dimensional VPFP system 

The CEM consists of inserting the following ansatz in 



-eVxF} , 



(71) 



up to terms of order e 2 . Since P — V x ■ E, E 
we may rewrite this equation as 



e 2 



(2tt)* 



P(x,t;e) + J2t n P in) (z,v;P). 



dP 
~dt 



J2e n F (n \P), 

n=0 

V = v + V x cb. 



(63) 

(64) 
(65) 

(66) 

(67) 
(68) 



Equations (£>2|) and (|33|) are obviously generalized to the 
three-dimensional case. The solvability condition for the 
first equation yields 

F (0) = V a -(V X ^P). 
Substituting F^ in we find 



E t + E(V X -E) = e {(V, • E) [{E ■ V X )E + E t ] 

+V X (V X -E)} + J, (72) 

where • J = 0, so that J is a solenoidal field. In the 
one-dimensional case, J = 0, E t — —EE X + O(e), and 
(f72|) becomes the Burgers equation d6j 



Then <fi obeys the Poisson equation 

= -P, 

and the normalization conditions are 
P{x,t)dx = 1, 

p( n \x,v;P) dxdv = 0. 



IV. CLOSURE OF EQUATIONS FOR THE 
MOMENTS 



Equations ( |6l| ) or (62) can also be obtained by ap- 
propriate closure of the hierarchy of equations for the 
moments of p, as we will show now. The equations for 
the first two moments of p, 



pdv, J — / vpdv 



e 2 



(27T)! 



Pv.V x V x . 



,(0) 



-V X P + PV X . 



i(0) 



P 



can be obtained directly from the VPFP system: 

Pt + Jx = 0, (73) 
-J -{<t> x P) x = cJ t + e— v 2 pdv. (74) 

If we take the x-derivative of the second equation and 
eliminate J by using the first one, we obtain 



5 



-ePt, 



dx 2 



v pdv. 



(75) 



Let us now close the hierarchy of equations by inserting 
the zeroth order ansatz (fhjl) 



v 2 pdv 



/2tt 



Pdv = (1 



iP 



(76) 



in (75). The result is the nonlinear wave equation |10| 



^-(^)*-{£r[(i 



%)P]-Pt> 



(77) 



to be solved together with (|4^) and ( }49| ) . For the electric 
field, £ = -<f> x , P = E x , and (|77j) yields 



E t + EE, r + eE tl 



e^[(l + E 2 )E x }=0, 



(78) 



after integrating once with respect to x and using the 
decay conditions to cancel the resulting constant. We 
can now show that this equation is compatible with our 
previously derived Burgers equation. Let us iterate (|7T 
to get 



Ett ~ —(EE x ) t — —E t E x — EE xt 
~ EE 2 + (EE X ) X E = (E 2 E X ) X . 



(79) 



Inserting this result in ([78]), we recover the Burgers equa- 
tion (|62[). In the limit as e — ► 0+, both the Burgers equa- 
tion and Equation ([78]) regularize shock waves and yield 
the same shock speed, (E + +E_)/2, for a shock connect- 
ing the field values E- and E+ However the inner 
structure of the shock is different in both cases and it 
would be interesting to study Equation ( j78|) for its own 
sake. 



V. CONCLUSIONS 

We have derived a reduced equation for the charge 
density (or the electric field) of the VPFP system in 
the high-field limit. The simplicity of the Fokker-Planck 
collision term allows us to explicitly find a zeroth order 
solution which is a displaced Maxwellian. This result 
is the basis of two classical singular perturbation proce- 
dures used to obtain the reduced equation: the Hilbert 
expansion and the Chapman-Enskog method. We find 
that the second method is ideally suited to provide a re- 
duced equation whose terms may be of different order. In 
the present case, we obtain a form of the Burgers equa- 
tion. On the other hand, closure of the equations for the 



moments of the probability density yields a nonlinear 
wave equation, which is compatible with the previous 
result, but may be of independent interest. Another 
interesting open problem is to use the CEM in the rigor- 
ous analysis of the high-field limit of the VPFP system. 
So far, rigorous analyses have used the less systematic 
method of closing moment equations. The trouble with 
using CEM seems to be that it is more "nonlinear" (al- 
beit systematic) than the other methods. 
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